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Abstract 

The next-to-minimal supersymmetric standard model (NMSSM) with an extended 
Higgs sector offers at least one Higgs boson as the Standard model (SM) like Higgs 
with a mass around 125 GeV. In this work, we revisit the mass spectrum and 
couplings of non-SM-like Higgs bosons taking into consideration most relevant 
constraints and identify the relevant parameter space. The discovery potential 
of these non-SM-like Higgs bosons, apart from their masses, is guided by their 
couplings with gauge bosons and fermions which are very much parameter space 
sensitive. We evaluate the rates of productions of these non-SM-like Higgs bosons 
at the LHC for a variety of decay channels in the allowed region of the parameter 
space. Although bb, tt modes appear to be the most promising decay channels, 
but for a substantial region of parameter space the two-photon decay mode has a 
remarkably large rate. In this study we emphasize that this diphoton mode can 
be exploited to find the non-SM-like Higgs bosons of the NMSSM and can also be 
a potential avenue to distinguish the NMSSM from the MSSM. In addition, we 
discuss briefly the various detectable signals of these non SM Higgs bosons at the 
LHC. 
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1 Introduction 


The recent discovery of the Higgs particle has completed the particle family in the Stan¬ 
dard model (SM), although more precise measurements are required to establish this claim 
strongly. On the other hand, there are various new physics models which can accommodate 
this Higgs candidate and are found to be consistent with the current measurements of Higgs 
couplings strength For example, the phenomenological models based on supersym¬ 

metric (SUSY) theory can interpret this Higgs candidate as their lightest Higgs boson [^[^. 
However, it is to be noted that to accommodate a Higgs boson of mass ~125 GeV in the 
theory of Minimal Supersymmetric Standard Model(MSSM), one requires a substantial con¬ 
tribution to the Higgs mass from higher order correction which is very close to the tree level 
value of it. Of course, this can be achieved in a very suitable region of SUSY parameter space, 
in particular with the contribution from the third generation squark sector along with some 
degrees of fine-tuning . 


The next-to-minimal supersymmetric standard model (NMSSM), an extension of MSSM 
with an additional Higgs singlet was proposed to solve the /r problem [^,11 . Interestingly, it 
has been found that this model can provide a SM-like Higgs boson of mass around 125 GeV 
more naturally (For example, see Refs. n2-20 and references therein), and as a result, after 


the discovery of the Higgs particle [^0, it has received a lot of attention. In this model the 
additional interaction term between the singlet and the doublet Higgs superfields contributes 
substantially to the Higgs mass, even at the tree level and thus reduces the required fine 
tuning significantly. The phenomenology of the NMSSM, more precisely of its Higgs sector, 
is quite rich due to the presence of this extra Higgs singlet and the corresponding interaction 
with the Higgs doublet. In the literature, there exist several studies carried out to study 
the NMSSM Higgs sector in the context of the LHG. The recent observation of the Higgs 
particle at the LHG and related measurements of its properties renewed the interest in the 
NMSSM leading to more focused and dedicated activities. The implications of the current 
observation of the Higgs boson at the LHG and the measurements of its properties have 
been extensively discussed in the context of the NMSSM by many authors. In the LHG 
experiment looking for the NMSSM Higgs signal in the current Run I data sets is one of the 
ongoing studies. In a very recent analysis, very light Higgs boson in the NMSSM was studied 
and from the non-observation of any signal, a bound on the Higgs production cross section 
times the corresponding Higgs branching ratio (BR) is presented 22 . 


The NMSSM contains seven Higgs bosons, three are GP even Higgs scalars (Hi, i = 1, 2, 3) 
and two are GP odd states, while there are two charged Higgs bosons. In this work we revisit 
the Higgs sector of the NMSSM by studying the interplay of model parameters and Higgs 
masses and couplings. We explore the possible implications in the NMSSM Higgs sector of 
recent Higgs discovery by scanning the model parameter space for a wide range requiring one 
of the GP even Higgs bosons SM-like Higgs with mass around 125 GeV. Performing a very 
comprehensive analysis taking into account various experimental and theoretical constraints, 
we delineate the region of parameter space which provides one of the GP even Higgs boson 
SM-like. Gorresponding to this allowed parameter space, we then discuss the phenomenology 
of, particularly, light neutral non-SM-like GP even and odd Higgs bosons at the LHG. It is 
observed that there exist a variety of decay channels of these light neutral non SM Higgs 
bosons with very reasonable BR depending on the parameters points. In this work, we study 
very systematically in detail the sensitivity of these BRs on model parameters and their impact 
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on light non-SM-like Higgs boson signals at the LHC. Moreover, in order to understand the 
feasibility of non-SM-like Higgs searches at the LHC, the rates of production of Higgs bosons 
in various decay channels are estimated for the center of mass energy t/s = 13 TeV, in 
the context of the interesting region of parameter space. In fact due to many parameters 
dependence of the cross sections and BRs of the NMSSM Higgs bosons, it is difficult to 
present systematically the variation of non SM Higgs boson production rates in various decay 
channels. Hence, we follow the strategy to present the ranges of these rates by computing the 
minimum and maximum values in various decay channels for a given Higgs boson mass and 
the parameter space. Naively, this estimation of rates provides hints about the detectability 
of non-SM-like Higgs bosons at the LHC although detail understanding of the background 
level is required to make a hnal conclusion. Interestingly, it is observed that the HR of non- 
SM-like Higgs in the di-photon channel is substantially large for a certain region of parameter 
space (l8,23 . Undoubtedly, this two photons channel appears as a striking feature of the 


NMSSM model. Perhaps, the di-photon channel can provide a robust mode to distinguish 
the NMSSM from the MSSM because of the fact that photon is a clean object to probe 
experimentally. In addition there are also other hadronic decay channels of non SM Higgs 
bosons which can be useful as well in looking for Higgs signal in the NMSSM. It is to be noted 
that the similar type of study is carried out, in particular in Ref. |24| where the main focus is 
to study only the pair production of Higgs. In this paper, we study the production of Higgs 
via all dominant channels focusing the range of masses for all Higgs bosons below 125 GeV, 
the interesting range in the context of the present LHC experiments. In order to discuss the 
Higgs signal we have very systematically discussed the HR of various decay channels of non 
SM-like Higgs bosons including the two interesting channels, gg and cc. More importantly, 
in this current study we obtain rates for various decay channels estimating the Higgs boson 
production cross section following the model based calculation implemented in SuSHi 
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The paper is organized as follows. In section 2 we discuss very briefly the relevant parts of 
the NMSSM Higgs sector. In section 3, describing inputs and constraints used while scanning 
parameters, we discuss the allowed ranges of the masses, the couplings and decay BRs of non- 
SM-like light Higgs bosons. Presenting the rates of the non-SM-like Higgs boson production 
in various channels in section 4, we discuss various Higgs signals in section 5. Finally we 
summarize our results in section 6. 


2 NMSSM Higgs sector 


As mentioned in the introduction, the NMSSM Higgs sector contai ns one extra gauge singlet 
Higgs superheld S in addition to two Higgs doublets (Hu,Hd) 
potential as. 
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making the super 


WnmSSM = WmSSM + A^HuHd + 


( 2 . 1 ) 


where Wmssm corresponds to the super potential in MSSM with only two Higgs doublets 
without a g term, A and k are dimensionless couplings. Recall that the primary goal in 
constructing the NMSSM is to generate the g term dynamically t o en sure its value around 
the EW scale, and thus solving the g problem. The last term in Eq. 2.1 with S^ is introduced 
to avoid Peccei-Quinn(PQ) symmetry 26 . Notice that this expansion respects explicit 
symmetry. In the MSSM, Wmssm depends on Hu and Hd, the Higgs super helds which couple 
with up-type quarks and down type quarks and charged leptons respectively. 
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Correspondingly, the soft-SUSY breaking potential for NMSSM Higgs sector is given by, 


Vsoft — 


|2 I 2 

I +“Hh 


|Hd|^ + mg|S|^ + 


AAASH^Hd + \kA^S^ + h.c 
o 


( 2 . 2 ) 


Here niHu, niH^, ms, Aa and are the soft breaking parameters. Notably, the parameters 
A, K, Aa are taken to be positive. The vacuum expectation value (VEV) of the singlet field, 
{S)=Vs, in turn generates the effective fj, term, //efr=Avs and then restricting it at the order 
of EW scale. The neutral components of other Higgs doublets, get VEVs as (Hu*^) = Vu and 
(Hd°) = Vd and the constraint m^ = ^(v^ + v^) makes one of the VEVs as a free parameter 
and it is parametrized as tan/3 = Vu/vd- 

At the tree level, the Higgs sector of NMSSM have 9 parameters i.e. 


A, K, tan/3, Aa, A^, m|. 


(2.3) 


But the minimization conditions of the scalar potential with respect to VEVs Vu,Vd and Vg 
reduce another three parameters, and hence, the NMSSM Higgs sector at the tree level is 
described by six independent parameters. 


A, K, tan,0, Hes, Aa, A«,, 


(2.4) 


unlike the MSSM where only two parameters mA, the mass of the pseudoscalar Higgs and 
tan/3 are required. More detail discussions and review on the NMSSM Higgs sector can be 
found in Refs. 


10-13,27,28 


The NMSSM Higgs sector contains 10 scalar degrees of freedom, out of which three are 
absorbed to give the masses of three gauge bosons, and Z leading to seven physical Higgs 
states. Expansion of Higgs fields around the three vevs v^, Vd and Vg yield various Higgs 
mixing terms in the Lagrangian of which real and imaginary parts constitute two 3x3 Higgs 
mass matrices corresponding to the CP-even and CP-odd scalars. The diagonalization of 
these mass matrices by orthogonal matrices provides masses of three physical Higgs bosons. 
Eor instance, the diagonalization of the 3x3 real mass matrix results in the masses of three 
CP even Higgs bosons (Hi,H 2 ,H 3 ) and the corresponding diagonalizing matrix relates the 
weak Higgs boson states, H" = (Hu^, Hd^, S^) with the physical ones as. 


Hi = SijHf; i,j = 1,2,3, 


(2.5) 


where Sy is the orthogonal matrix diagonalizing the 3x3 CP even Higgs mass matrix. Here the 
CP even mass eigenstates (Hi) arranged in increasing order of masses i.e muj < mHj < mHg. 
Consequently, the physical Higgs boson states turn out to be admixtures of both doublet 
(Hu, Hd) and singlet (S) states and naturally the corresponding mixing angles affect their 
couplings to fermions and gauge bosons. Similarly, the imaginary parts of Higgs fields form a 
3x3 mass matrix in the basis A" = (Hf,, H)^, S^) and is being diagonalized by an orthogonal 
matrix as. 


Ai = PijAf, (2.6) 

Ai stands for three CP odd states, (Ai, A 2 , G°), where is the Goldstone boson. Defining a 
suitable 2x2 matrix, the Goldstone boson state G^ can be rotated away leading only to two 
physical states Ai and A 2 which are admixtures of doublet and singlet states determined by 
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the components of the orthogonal matrix Py. The mixing angles i.e the matrix elements of 
the orthogonal matrices depend on the model parameters given by Eq ]2.4| [^,27 28 and they 
play an important role in determining the couplings of the Higgs bosons with fermions and 
gauge bosons making those very much parameter space sensitive. For example, the couplings 
of CP even Higgs bosons Hj to fermions, of both top and bottom type, are given by p p7|28l 
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and similarly for CP odd scalars. A,, 
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( 2 . 8 ) 


Here Su, Si 2 , (Pii, Pi 2 ) determine the amount of doublet composition in physical CP even(odd) 
Higgs boson states. Thus due to the presence of mixing angles, Higgs couplings either get 
suppressed or enhanced. For instance, if any of the Higgs state (Hi, Ai) be dominated by singlet 
components i.e Sii,Si 2 ~ 0, Pii,Pi 2 ~ 0 then the corresponding Higgs boson couplings with 
fermions and gauge bosons are severely suppressed. The ratio between the Higgs couplings 


given in Eq. |2.7| |2.8| and its SM counter part is just a scale factor, called reduced coupling 
depending on mixing angle and tan/3, e.g Cu = Su/sin/3 is the reduced coupling of Hi with 
up-type quark. Moreover, the orthogonality property of the mixing matrices(S,P) predicts a 
sharing of couplings among the corresponding Higgs boson states. For example if any of the 
CP even Higgs boson be SM like i.e Su or Si 2 ~ 1, then the corresponding couplings of other 
CP even Higgs bosons with fermions are heavily suppressed. The coupling of CP even Higgs 
bosons with gauge bosons (V=W,Z) are given by 


13,28 


gHiVV = gH^vv X 6; 6 = cos/3Si2 + sin/3Sii, 


(2.9) 


and because of orthogonality property of matrix S, one can show that the mixing angle follows 
the relation, 




( 2 . 10 ) 


As stated earlier this sum rule implies that if any one of the Higgs boson couples dominantly 
with a gauge boson, then the same couplings of other Higgs bosons are severely disfavoured. 
The mixing angle is constrained from the non observation of NMSSM Higgs bosons in LEP 
experiment |29| . 

In NMSSM the upper bound of lightest Higgs boson mass at tree level is well known to 


be as 10,11,27,30,31 


m^, < M| cos^ 2/3 + A^v^ sin^ 2/3, 


( 2 . 11 ) 
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where the first term is similar to the upper bound (=Mz cos 2/3) in MSSM and the second term 


originates due to the interaction between singlet and doublet Higgs field via A term, in Eq. 2.1 


Remarkably, the contribution of the second term pushes this upper bound substantially much 
above Mz cos2/3, of course, depending on the values of A and tan/3. Clearly, high (low) value 
of A (tan/3) is preferred to obtain a large value of the lightest Higgs mass at the tree level. 
Hence, this tree level NMSSM contribution to lightest Higgs mass enable to accommodate 
easily a SM-like Higgs mass of the range ~ 125 GeV. 

As pointed out earlier, the Higgs sector also receives higher order perturbative corrections 
contributed by squark masses and trilinear A-terms, and as well as by the parameters of the 
gaugino and Higgsino sectors. A very tiny contribution may come from slepton sectors also. 
The total extra contribution due to these corrections lead further enhancement of Higgs mass. 
The calculations of these higher order corrections exist in the literature C>(at«s + afaCts) (3^ 


41 and the NMSSM-specific two loop corrections can also be found in the Ref. [^. Including 
this additional correction to the mass, this upper limit of lightest Higgs boson mass can go 
up to ~ 140 GeV with the constraint that A to be within the perturbative limit [42[|43| . 

Consequently, the higher order correction to Higgs mass relates the other sectors, in 
particular third generation of squarks. Therefore to predict Higgs mass, one needs to take 
into consideration of involved model parameters, including very sensitive third generation 
sector. On the other hand, perhaps, the well determinations of Higgs sector predict about 
the range of masses and mixing angles of third generation squarks, namely top squarks [6,44 


However, it is to be noted here that the different spectrum generators for the NMSSM predict 
different Higgs boson masses separated by few GeV for a given set of input parameters and 
schemes 1451. 


As described before, in NMSSM the Higgs sector has a non trivial dependence on a 
set of parameters Eq|2A| [28] which eventually affect the NMSSM Higgs phenomenology. 
Recall that generating /r term dynamically leads /ieff ~ Avg, which also is connected with 
chargino/neutralino masses. Hence, the LEP limit on light chargino mass ( > lOOGeV) 
[46] predicts a bound, A > //eff/vg [^. Apparently, this bound restricts A not to be very 
small unless vev Vg becomes too large, which is not also phenomenologically favourable as the 
other vevs, Vu,Vd are expected to be also around the EW scale. However, there are certain 
variations of the NMSSM model, for instance, the semi-constrained NMSSM where soft masses 
of scalars, gauginos and trilinear A-terms for third generation fermions are assumed to be the 
same at the grand unified scale [^. In this type of scenario, the values of the A and k can be 
very small, even much below of the 0(0.1) which is completely compatible with the existing 
constraints which will be discussed later. On the other hand, running of A and k from high 
scale of Grand Unified Theory(GUT) to EW scale through renormalization group equation 
and the requirement of perturbative nature of couplings, a bound on them can be derived 


as, 


A^ + < 0.6. 


( 2 . 12 ) 


3 Scanning of NMSSM parameters 


In this section, we try to identify the region that offers one of the CP-even Higgs boson as 
the SM-like (Hsm)- by carrying out a comprehensive random scan of NMSSM parameter 


space. We use the package NMSSMTool4.3.0 48,49 which calculates for a given set of input 
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parameters, the masses, the couplings and the decay branching ratios of all the Higgs bosons 
along with the spectrum of other SUSY particles. This code also systematically checks the 
consistencies of the parameters against various theoretical as well experimental constraints. 
Experimental constraints include measurements of various low and high energy observables 
predicted within the framework of the NMSSM. 

In our scan, we set the following ranges for the relevant parameters. 


0.1 < A < 0.7; 0.1 < K < 0.7; 0 < < 2 TeV, -2 TeV < < 200 GeV; 

1.5 < tan/3 < 30; 100 GeV < /Tes < 2 TeV. (3.13) 


Here the choice of parameters are phenomenologically motivated which expected to provide 
a very different scenario than the MSSM and testable at the future LHC experiments. The 
range of parameters are chosen in such a way that it provides one of the CP even Higgs boson 
to be SM-like with a mass around 125 GeV and other singlet like Higgs bosons are not too 
massive. In our scan very small values of A are avoided to restrict //g// and Vg at the level 
of electroweak scale. However, as mentioned in sec.2, the values of A and k can be very 
small even much lower than our considered ranges, in the context of the certain variation of 
the NMSSM [47] . In view of this, our chosen scenario of the paratmeter choices is not very 
general one, represening only a subset of the parameter space. The g-2 constraint favours 
the positive values of ^ieff■ However, this in not a very general set up of parameters while 
more wider ranges of parameters are also accessible in the context of certain variation of 


NMSSM 47,50,51 


The soft mass parameters for the left and the right handed squarks from the first two 
generations, to which the Higgs boson mass are not so sensitive, are set to, 


MQi /2 = ^Ui /2 = Mdi/ 2/3 = ITeV, (3.14) 

while for the third generation squarks, known to have a large effect on the Higgs mass, these 
are varied over the range. 


Mq 3 = Mu 3 = 300 - 3000 GeV. (3.15) 

to consider a wide range of values from them. However, as we know, because of radiative 
contributions to the mass of the Higgs boson, the masses of the top squarks and the corre¬ 
sponding mixing are important and hence A-term plays an important role. With a goal to 
achieve a GP even Higgs state close to the mass of the SM Higgs, we vary |At| from small to 
large values as. 


At = -4 TeV to + 4 TeV, (3.16) 

and setting other third generation trilinear A-terms as, 

Ab = 2 TeV; and Ae 3 = 1 TeV, (3.17) 

while keeping all other remaining trilinear A terms to zero. While scanning parameters, we 
test all constraints including theoretical ones pertaining to vacuum stability. Moreover, the 
precision measurements at the LEP experiments restrict the new physics models through the 
measurement of Mw, sin^ 0eff) the p parameter and also via the invisible decay width of the Z 
boson. In addition, the constraints due to various low energy observables from flavor physics 
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such as B —)• Xs 7 , Bg —)• and B —>■ and also the mass differences AMj and AMg are 

checked in the process of scanning of the parameters. Various existing bounds on sparticle 
masses and production cross sections that are obtained from LEP, the Tevatron and from 
the LHC experiments are also imposed. The details of all these constraints can be found 
Notably, among various experimental constraints, the most crucial ones are from the 
5^ data and the muon anomalous magnetic moment |54|. In this model the lightest 


m 


52 


Planck 

neutralino(x5) is assumed to be the lightest-supersymmetric particle (LSP) and is also a DM 
candidate. NMSSMTools interfaced with micrOMEGAs (55,56 , calculates the relic density. 


Recent data from the PLANCK experiment |53| concludes the relic density to be, 


D^h = 0.1187 ±0.0017. 


(3.18) 


In the NMSSM, the DM solution is tested with the measured data as shown above at the 
10% level. The most natural solution to the DM problem is via the bulk annihilation channel 
which requires lighter electroweak gauginos and sleptons. The impact of DM constraints on 
the NMSSM are discussed by many authors [16,47,56-66 . The parameters in the gaugino 


sector. Ml and M 2 , the U(l) and the SU(2) gaugino masses, respectively are expected not to 
have any large effects on Higgs masses, but are important in calculating the dark mater relic 


density, Eq. 3.18 Therefore, we vary the values of gaugino masses Mi, M 2 while keeping M 3 , 
which is very close to the value of the gluino mass to a fixed value. 


50 GeV < Ml < 1 TeV; 50 GeV < M 2 < 1 TeV; M 3 = 1.2 TeV. 
The precisely measured value of the muon anomalous magnetic moment, 

a -feciL 

O'// — ^ 


(3.19) 


(3.20) 


is considered to be a very strong constraint for new physics. The total SM contribution to a^ 
is estimated with an uncertainty predominantly due to the hadronic contributions 67 . The 


discrepancy between the measured value and the one predicted by the SM is found to be 54 , 


<5a^ = (28.7 ±8.0) x 10“^°. 


(3.21) 


This implies the measured value is 3 cj away from the SM predicted value. In SUSY model, 
there are additional contributions due to the presence of sparticles in loops which provide 
an explanation for this excess 68-71 . In the NMSSM, the effect of g^-2 constraint on its 
model parameter space has been discussed in (7^. As we know, the main SUSY contribution 


to a^ comes at the one loop level involving smuon-neutralino and sneutrino(muon)-chargino 
diagrams. However, this SUSY contribution at the one loop level is determined by the sign 
of the /i-term and the value of tan (3. In order to have a reasonable contribution to da^ from 
these additional one loop SUSY diagrams, the masses of the sparticles, in particular of the 
smuons are favoured to be as light as ~ 100 GeV and tan/3 ~ 10-15 [^. The slepton masses 
are not directly coupled with the Higgs sector, but plays a crucial role in determining the 
anomalous magnetic moment of the muon as discussed before. In order to have a parameter 
space consistent with 5^-2 constraint we assume low values of the first two generations of 
sleptons 72 i.e. 


Mlj^ 2 = 100 GeV; 


Mj 5 i 2 — 100 GeV. 


(3.22) 
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Lower 

Upper 

Ct 

0.97 

2.28 

Cb 

0.00 

1.23 

Cv 

0.66 

1.23 

Cg 

0.52 

1.07 

C-y 

0.67 

1.33 


Table 1: The allowed ranges of reduced couplings of the SM-like Higgs with fermions and the 
gauge bosons and also the effective couplings with the gluons and the photons at 95%C.L. 
These are obtained by CMS experiments [^. 



In our scan, we take into account the impact of g^-2 constraint which has some effects on the 
parameter space as will be discussed later. 

In the NMSSM, out of three CP even Higgs bosons, one is required to be SM-like. In 
addition, the couplings of the SM-like Higgs boson to fermions and the gauge bosons should be 
consistent with the current measurements by both the CMS and the ATLAS experiments . 
As shown in EqU^ |2.8| and |2.9[ the couplings of the Higgs boson with fermions and the 
gauge boson are basically a scaling of the corresponding SM couplings by a factor called the 
reduced couplings which are essentially the ratios of gHiff/gH^(= Cf) or gHiVv/gH^v(~ Ey)- 
Similarly, one can also have reduced effective couplings of Higgs with gluons and photon, 
say Cg and C-^, respectively. At the LHC, both the ATLAS and the CMS collaborations 
examined the the Higgs couplings following the strategy discussed in Ref. |75| . By performing 
a very detailed analysis, the allowed ranges of the reduced couplings are reported at 95% C.L. 
In Table allowed ranges of these reduced couplings measured by the CMS experiment 
are presented. Note that the CMS and ATLAS experiments do not measure the Higgs 
couplings directly. However, by measuring the signal strength, the factors by which the 
SM couplings are scaled can be constrained. Here these scale factors are defined to be the 
reduced couplings which are constrained by both the LHC experiments. It is to be noted that 
the SM Higgs production cross sections in various related modes are used in deriving these 
limits. In our study, while scanning the parameter space, we impose these constraints on the 
reduced couplings for the SM-like CP even Higgs boson. Similarly, from negative searches of 
Higgs at the LEP experiment, the reduced couplings of the SM-like Higgs boson with W/Z 
(Cv = C) a-lso get constrained 29 and our scan has been subjected to these constraints as well. 
Furthermore, in the NMSSM, possibly the SM-like Higgs boson can decay via non-standard 
modes which will be discussed in later section. The total BR in these non SM modes, namely 
BRbsm of the SM like Higgs boson is also restricted from above from the current Higgs data. 
The CMS measurements set the upper limit of non-SM branching ratio of the SM-like Higgs 
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to be [^, 


BRbsm < 0.57 at 95%.C.L. 


(3.23) 


This limit is incorporated in our parameter scan. 

Eventually, after playing with the NMSSM model parameters, we identify the region of 
parameter space allowed by all theoretical and experimental constraints as described.above, 
which yields one of the relatives light CP even Higgs bosons, either Hi or H 2 , to be SM-like 
with its mass in the range, 

mHi/mHa ~ ™Hsm = 125.02 ± 3 GeV. (3.24) 

It is always observed that the heaviest CP even neutral Higgs, H 3 tends to be massive ~ 200 
GeV or more. Thus, we focus our attention on two distinct regions presenting the following 
two scenarios: 

• Gase A: the second lightest CP even Higgs, H 2 ~ Hsm- 

• Case B: the lightest CP even state. Hi ~ Hsm- 

In the subsequent sections we discuss various features of these two cases and the relevant 
phenomenology in the context of ongoing LHC experiment with a center of mass energy of 
13 TeV. 

3.1 Ccis0 A: H 2 ~ Hsm 

In this section we discuss the salient features of the NMSSM Higgs sector corresponding to 
the parameter space which leads to the second lightest CP even Higgs boson, H 2 , as SM-like. 
Here, the lightest CP even Higgs(Hi) lighter than H 2 by definition. In addition the lightest 


g-2 g-2-^DM-^LHC ■ 



Figure 1: Range of niA^ and mH^ for case(A), allowed by all theoretical and experimental 
constraints as discussed in Sec. 3 plus — 2 constraint (cyan) and then adding the LHC and 
DM constraints (blue). 


CP odd Higgs boson(Ai) is also lighter than H 2 in a certain region of parameter space. Both 
Hi and Ai are dominantly singlet like. In Fig. we present the possible ranges of mH^ and 
lUAi when H 2 is SM-like 24 . In this figure the allowed points (cyan) are result from imposing 
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only the gij_-2 limit (Eq j3.2l] ) along with all other theoretical and experimental constraints as 
described in the previous section. Subsequently, adding constraints from the DM relic density 
(Eq. 3.18) and from the measurements of the Higgs couplings at the LHC (Table[^, the more 
restricted region(blue) is obtained. It is interesting to note that, even after including all such 
constraints. Hi or Ai may turn out to be very light, ~ 20 GeV, although there are only few 
points [76 77 . Moreover, this figure clearly indicates that for a good fraction of the parameter 


space, both mH^ and mA^ are found to be less than the half of the mass of H 2 , thus allowing 
the latter decay in the following ways. 


Ho 


HiHi, 

AiAi. 


(3.25) 


It is to be noted that these lighter singlet states Ai and Hi are not yet ruled out by any collider 
experiments. These singlet like states escape detection in collider experiments because of their 
suppressed couplings to gauge bosons and fermions. It can be realized by looking at Figsj^and 
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Figure 2: Reduced couplings of Hi with gauge bosons(V=W/Z) and fermions (up and down 
types) for case(A). Values are allowed by mentioned constraints ( conventions are the same 
as in Fig. j^) 

which present for a wide range of tuhi i reduced couplings of Hi and Ai, respectively, 
with gauge bosons, up and down-type quarks. At lower masses for Hi, allowed values of 
these reduced couplings initially appear to be very small and then go up with increasing mH^ • 
However, LHC constraints, in particular restrictions on reduced couplings of the SM-like Higgs 
boson as shown in Tablej^ favor smaller values, except for the case which receives some 

enhancement for higher values of tan/3 (see Eq. 2.7). In case of Ai, the pattern of variation 
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of reduced the couplings with uiai is a little different since Ai physical state contains a finite 
fraction of doublet component and the reduced couplings increase for higher values of tan /3. 
As observed in Fig.[^ the reduced coupling of Ai to d-type quark is somewhat larger than its 
coupling to u-type quark. Consequently, this behaviour of the reduced couplings significantly 
affects the Higgs phenomenology in the colliders. These are discussed in the next section. 
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Figure 3: Reduced couplings of Ai with quarks (up and down types). Other conventions are 
same as in Fig. 

3.2 Case B: Hi ~ Hsm 

Following a similar strategy as described above, we obtain a substantial region of parameter 
space with the lightest CP even state Hi as the SM-like Higgs boson. Obviously, in this 
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Figure 4: Allowed range of and mH 2 for case (B). Constraints are the same as in Fig. 

scenario, mHj ^ niHi ~ aiHg^ and the mass of Ai can be lighter or heavier than Hi. In fact, 
a signihcant region of parameter space exists even after imposing all constraints where it is 
observed that Ai is light, around ~ 50 — 100 GeV and it can turn out to be even lighter 
~ 10 GeV for few points( see Fig. [^, where allowed ranges of ruAj and mH 2 are presented 
for this scenario. This figure indicates that the available range of mA^ can be lighter than 
m/i^/2, whereas may extend to very large (~ 500 GeV) values. Hence the following 
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decay channels, 


H2 ^ HiHi,AiAi, 

Hi ^ AiAi, (3.26) 

open up with a reasonable branching ratio depending on the region of the parameter space. 
These decay channels may potentially provide additional indirect sources of productions of 
Hi as well as of Ai at the LHC mediated by triple Higgs bosons couplings. 
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Figure 5: Reduced couplings of H 2 with gauge bosons (V=W/Z) and quarks( up and down 
types) for case(B). Constraints are the same as Fig. 
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Figure 6: Same as in Fig. but for case(B). 









As before, in order to understand the phenomenology in this scenario of non SM singlet¬ 
like H 2 and Ai, we investigate the coupling strengths of these Higgs states. Since Hi is 
SM-like, and it predominantly couples to gauge bosons, obviously the sum rule (Eq. |2.10[ ) 
predicts suppression of the corresponding reduced couplings of H 2 , as shown in Fig. [^(bottom 
panel). Notice that the reduced couplings of H 2 with d-type quark are reasonably large, ~2- 
10 depending on the masses, but additional DM and LHC constraints restrict those to lower 
values. The variation of the reduced couplings of Ai is also studied and is shown in Fig. 
The reduced couplings of Ai to the u-type quark is almost negligible whereas with the d-type 
quark it has finite values. Moreover, because of the absence of Ai coupling to gauge bosons, 
direct production of Ai in lepton colliders and via vector boson fusion in hadron colliders are 
suppressed. Hence, the low masses of Ai are not excluded by LEP experiment. However, in 
hadron colliders Ai can be produced via gluon gluon fusion or b -b annihilation and due to 
the suppressed couplings of Ai with u-type quarks, the latter process is the comparatively 
dominant one. 

Armed with this knowledge about the spectrum of allowed masses and the couplings of 
the light neutral non SM Higgs bosons for the two scenarios, (A) and (B), we now discuss 
their decay patterns and followed by their phenomenological implications at the LHC. 

3.3 Higgs decays: Hi,Ai,H 2 

As pointed out earlier, the NMSSM Higgs bosons offer a rich phenomenology at colliders owing 
to their very diverse decays in various channels including some non SM modes. Furthermore, 
the BR of each Higgs decay modes and hence, the sensitivity of the corresponding signal, are 
very much parameter space dependent because of the presence of complicated admixtures of 
physical Higgs states making them either doublet- or singlet-like. In this work, in order to 
predict NMSSM Higgs signal at the LHC in the context of the scenarios presented in case(A) 
and case(B), we revisit various Higgs decay channels by studying correlations among the BRs 
and the parameter space. We also focus on a few interesting Higgs decay modes which can 
lead to robust signals of the NMSSM Higgs bosons, in particular, for the lightest CP odd 
and the CP even Higgs bosons. In the present context, we discuss features of BRs of various 
decay channels only for two CP even (Hi,H 2 ) and the lightest CP odd (Ai) Higgs bosons. 

• Hi decays : As mentioned before, for certain regions of the parameter space in the scenario 
represented by case(A) with H 2 as the SM-like Higgs boson. Hi and Ai can be very light and 
can have masses below 50 GeV, see Fig[2 The decay Hi —)• AiAi will be the dominant one 
whenever it is kinematically allowed. The corresponding decay width is primarily determined 
by triple Higgs Hi — Ai — Ai coupling mediated by singlet components and approximately 
proportional to KVg, thus making the decay width very large in comparison to other fermionic 
and bosonic channels. It is found that for certain set of parameters, the BR (Hi —)• AiAi) 
goes to 90% or more. As expected, this triple Higgs boson coupling receives corrections [78|7^ 
which are not taken into account in this present calculation. There exists a subset of parameter 
space where Hi decays dominantly to bb channel with BR close to 90% and rest to the rr 
channel (~ 10%) as shown in Fig. j^left) for various masses of Hi. Although the Hi state 
is dominated by the singlet component, presence of a finite fraction of doublet component 
(Si 2 , composition of Hj) leads to a sizeable coupling of Hi to d-type fermions (Eq |2.7[ ). As 
a consequence, the partial widths of Hi —)■ bb, rr modes turn out to be dominant leading to 
enhancements of respective BRs. 
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Figure 7: BR of Hi in different channels for case(A) and allowed by all constraints including 
g^-2, LHC and DM. 


large 80 


BR is found for the 77 channel 81 


On the other hand, for certain region of the parameter space (where > 60 GeV), 
the coupling Chiuu becomes finite (see Fig.2) resulting in an enhancement of the decay width 
of Hi —)• gg (via top quark), cc channels leading a little suppression of bb decay BR as 
demonstrated in Fig. (right), provided Higgs to Higgs decays are kinematically forbidden. 
Notice that for the mass range around mui ~ 60 GeV or above, the gg, cc channels are quite 
A finite fraction of about 4-6% of Hi decays to WW for higher mui and also ~2-3% 

Note that the pattern of BR as shown in Fig. [^left and 
right) correspond to two different regions of the parameters space. For case(B), Hi is SM-like 
and decays via various SM channels with almost the same rate as in the SM. In addition, 
it can also decay to some non SM modes, such as Hi —?■ AiAijX^x)* whenever kinematically 
allowed. Note that, however, as mentioned earlier, the total BR to these non SM channels 
of SM-like Hi is restricted by an upper bound, as given by Eq. 3.23 Perhaps, looking for 
these non SM decay channels of the SM-like Higgs can be a potential avenue to confirm the 
existence of this type of SUSY scenario. 

• Ai decays : For the scenario presented by case (A), the lightest GP odd Higgs boson Ai 
primarily decays to fermions due to the absence of tree level couplings with gauge bosons. As 
shown in Figj^ for lower mass range, it decays dominantly to bb with BR ~90% and then to 
sub-dominant rr channels. However, if decay modes such as, Ai —)■ X 1 X 1 JX 1 X 2 accessible 
for a heavier Ai ( > 150 GeV), then those turn out to be the dominant ones due to the 
enhancement of Ai couplings with singlino-Higgsino components in y)*, X 2 states leading to 
larger width. Here, x? and X 2 are the lightest and the second lightest neutralinos respectively. 
In this scenario, the BR of Ai in the di-muon final state is very tiny (~ 10““^) as shown in 
the same figure along with other sub dominant decay modes such as, Ai —)• 77 , gg, Z 7 , ZHi. 
For the scenario where Hi is SM-like i.e.in case(B), the BRs of Ai to various final states are 
depicted in Fig. Interestingly, note that BR(Ai —>■ 77 ) can indeed be very large (~100%) 
and even greater than the BR for bb channel, in certain regions of the parameter space and 
for a wide range of myi^(left). This interesting decay pattern of Ai in the di-photon channel 
can be understood by more careful investigation of the structures of Ai couplings to various 
particles. 
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Figure 9: Same as in Fig. [^but for case(B). 


Here, note that Ai is dominantly singlet like resulting in all of its available fermionic decay 
channels suppressed. Primarily the off-diagonal elements in the Higgs mass matrix determine 
the composition of Higgs boson states. For example, in case of Ai when the off-diagonal term 
of pseudoscaler Higgs mass matrix vanishes for a certain combination of the parameters, then 
Ai becomes completely singlet like. From the 2x2 Higgs CP odd mass matrix, one can 
conclude that the vanishing of off-diagonal term leads Ai purely singlet dominated which can 
be translated to the relation such as Ax ~ it occurs for a wide region of parameters. 

Moreover, Ai does not have any tree level coupling with the vector bosons. In this kind of a 
scenario, Ai can still couple to photons at one loop with charginos in the loop. The Higgsino 
like charginos favourably couple to a singlet like Ai state and enhance the partial width of 
Ai —)• 77 mode. Therefore, for the regions of the parameter space where both the charginos 
are Higgsino like, BR(Ai —)• 77 ) turns out to be very large for a wide range of mA^ as clearly 

Evidently, this di-photon final state appears to be one of the striking 
features of this model and can be exploited not only to discover Ai at the LHC, but also 
to discriminate the MSSM from the NMSSM. It will be discussed later in more detail. In 
addition, Ai —)• rr mode is also present with a reasonable BR(~ 10%) which can also provide 
a clean signal of Ai. The other sub dominant decay modes of Ai are presented in Fig. 
(right) where the dominant one is the invisible decay channel, Ai —?• XiXi- The BR of Ai in 
the muonic channel, Ai —?■ /r/r is of the same order (~ 10 “^) as before along with other sub¬ 
dominant channels like Ai —>■ gg, Z 7 . Note that, both the plots in Figs, [^and [^presenting 
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the BRs of Ai and Hi respectively correspond to two different regions of the parameter space 
for each cases. 

• H 2 decays: In the scenario described in case(B), H 2 can decay to a pair of SM-like Hi as 
well as to a singlet like Ai pair i.e. H 2 —)• HiHi, AiAi. If kinematically allowed, these are 
expected to be dominant(~ 80-100%) is clearly shown in Fig. 10 (left). Since, Hi is SM-like, 
the decay channel, H 2 —HiHi is likely to contribute to the rate of SM di-Higgs production 
at the LHC. Similarly, in Fig|^ (right) we present BRs of various decay channels of H 2 for 
case(A). 

In this case, notice that for a good fraction of parameter space the mass of H 2 Hsm is 
heavier than that of Hi and Ai(see Fig 1), and BR(H 2 —)• HiHi) is found to vary between a 
few percent to ~40% where as BR{H 2 —)■ AiAi) is enhanced to 30%. Note that, as before, 
these are the non SM decay modes of the SM-like Higgs H 2 and also allowed by the constraint 
in Eq. 3.23 However, the decay channel H 2 —)• HiHi is expected to contribute to the visible 


signal of the SM-like Higgs boson production. Hence, it is worth examining the contribution 
of this decay channel in the present Higgs data collected at LHC Run I experiments, which 
eventually may constrain the model |83| . Notice that in this scenario, case(A), H 2 has a 
reasonable BR (~10-90%) to a pair of lightest neutralinos contributing to its invisible decay 
width 84 . In addition, H 2 has many other sub dominant decay modes which can provide 


interesting signals, in particular for the scenario (B). For example, decays to SUSY particles 
such as a pair of neutralinos, H 2 —)• X 2 X 2 ) X 2 X 3 ) provide clean signals because of leptonic 
decays of the neutralinos. We found that in this regions of the parameter space, H 2 —?• ZAi 
is very tiny. 
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Figure 10: BR of H 2 for case(B)(left) and case(A)(right) subject to all constraints as in 


Fig0 


4 Higgs production at the LHC 

In general, the NMSSM Higgs bosons can be produced via conventional mechanisms, i.e by 
the dominant gluon-gluon fusion (ggF), whereas the next sub dominant mode is the associated 
productions with a pair of b quarks. However for a very large value of tan/3, the associated 
production mode appears to be the dominant one. We reiterate here that in the NMSSM, the 
Higgs couplings to fermions and the gauge bosons are essentially the SM couplings scaled by 
the respective reduced coupling factors which can vary widely from small (<C 1) to moderately 
large values(~ 0.5). It implies that the Higgs boson production cross-section in the NMSSM 
can be obtained from the corresponding SM cross-section by appropriately rescaling them by 
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the respective factors. Hence the following naive strategy can be adopted to estimate the 


NMSSM Higgs boson production cross section as 24,85 


crggF(</>) = crl^i(p) X C' 


gg’ 


(4.27) 


where (j) stands for any of the CP even or odd Higgs boson states. Here crggp{(p) is the NMSSM 
Higgs boson production cross section via gluon-gluon fusion and o'g^((/>) is the corresponding 
SM cross section for the same Higgs boson mass. Cgg is the effective reduced coupling gg(f> via 
loops comprising heavier quarks and squarks. The SM Higgs boson production cross-section 


through ggF is known at the level of next-to-next-to leading order(NNLO) QCD 86-92 . We 


compute these SM Higgs boson production cross section via ggF using SusHi-vl.5.0 93 which 


takes into account the QCD contribution up to NNLO mediated by only the SM particles. 

However, it is to be noted that in SUSY model, the ggF channel may receive a sizeable 
contribution from loops involving the corresponding superpartners, in particular, the third 


generation squarks, the stops and the sbottoms 94 . Recently, in the context of the NMSSM, 


the neutral Higgs boson production cross-section via ggF is computed including NLO contri¬ 
butions of the squarks and the gluino in addition to the electroweak corrections involving light 
quarks and it is implemented in SusHi 25 . These calculations are based on the expansion 
of terms of inverse of heavy SUSY particle masses 95- 97 . Obviously, numerical estimation 


using this calculation is expected to be more reliable than the approximate method described 
Hence, in this work we use this calculation implemented in SusHi-vl.5.0 


in Eq. 4.27 


93 


to compute the neutral Higgs boson production cross sections in the NMSSM. Thus, the 
production rate of Higgs decaying to any of the pair of state X can be estimated as. 


Rxx = o-ggF((/>) X BR(()) —> XX), (4.28) 

where (j) —)■ XX collectively represents the Higgs decay channels with XX = bb, rr, 77 , gg, cc, /r/i, 
WW, ZZ etc. For a given decay channel and a set of parameter values, BR((/) —)■ XX) is ob¬ 
tained from the NMSSMTools. aggF{(j)) is computed directly using the NLO calculation 
in the NMSSM implemented in SusHi. On the other hand, for the sake of comparison, we also 
estimate the cross section for the same Higgs boson production with the same mass following 
the approximate method as described in Eq. |4.27[ R is found that following this prescrip¬ 
tion the calculated CP even Higgs boson production cross sections are 30-40% higher than 
the value obtained using a more exact calculation of SusHi |25| . This can be attributed to 
the fact that in the NMSSM, the production cross sections are estimated at NLO including 
electroweak corrections, where as in Eq. |4.27| NNLO QCD result without any SUSY loops are 
used along with Cgg estimated at NLO. 

The cross-section for the next-dominant production mode, can be obtained by exclusive 
Higgs boson production in association with a pair of b quarks , 

gg —)• hh(j) (4.29) 


25 


However, this cross section can be obtained also by computing bb —)• ()> which is a good 
approximation within the higher order correction (^|100 . Note that the SusHi calculates this 
bb —)■(/> at NNLO QCD. The rate Rxx corresponding to this production channel is obtained 
by multiplying the effective reduced hhcj) coupling factor with the SM production cross section 
obtained from SuSHi 93 and the respective BR((/> XX). In our calculation we use the 
CTIO (M for parton distributions function and both the QCD scales are set to for ggF, 
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whereas the factorization scale is fixed to 0.25M^ only for bb —)• (?!) production. Needless to 
say that the rate of the Higgs production in any of the decay channels is very much sensitive 
to the model parameters owing to the presence of mixing angles in the Higgs couplings, see 
Eq. |2.7| and |2.8[ 

Indirectly, Hi and Ai can also be indirectly produced in pairs through H 2 production via 
ggF and its subsequent decays provided mHi,mAi < mH 2 / 2 . The corresponding production 
rate can be obtained as. 


a. 


ggF 


(H 2 ) = aggF(H2) X BR(H 2 ^ </-(/)) 


(4.30) 


where ()) = Hi, Ai. As before, the crggF(H 2 ) is computed from NLO calculation in the NMSSM 
framework using SusHi |25|. 


Now we present the light neutral non-SM Higgs bosons production rates Rxx (Eq. 4.28) 
via both the production modes, ggF and bb —)• (/> for various decay modes corresponding to 
the regions of parameter space presented by case(A) and case(B). As pointed out earlier, the 
opening up of many decay modes of the non-SM-like Higgs bosons offers a variety of signals 
with a widely ranging rates. In order to demonstrate, we compute the minimum and the 
maximum values of those rates, Rxx for those region of the parameter space for a given mass 
of the Higgs boson and the decay channel. Naively, such regions of the parameter space which 
leads to the minimum and maximum cross sections for a given Higgs mass and decay channel 
are identified by the corresponding values of the product, Cgg x BR(0 —)■ XX). We assume 
that the parameter region which gives the maximum (minimum) value of this product yield 
a maximum (minimum) value of the cross section for the given Higgs boson mass and the 
decay channel. In Tables we present minimum and maximum rates for non-SM-like 

Higgs bosons in various decay channels for the scenarios case(A) and case(B). For example, 
in Table the results are shown for the non-SM-like Hi for case(A). Rates are presented for 
the Higgs production via both ggF {Rxx{gg)) and associated production, bb —)• 4>{Rxx{bb))- 
In this table, the first row shows the selected values of mui extending up to 120 GeV, while 
the first column indicates the rates (Rxx) for Hie decay channel XX of which minimum and 
maximum values are presented. For instance, with a choice of lower value of mH^ = 13 GeV, 
the maximum rate of Hi production in bb channel via ggF is about 145x10^ fb, where as the 
same for 77 and rr channels are about 1.4 fb and 14 xlO^fb, respectively. As seen in this 
Table that the rates (Rxx (bb)) via associated production Higgs are also not negligible. Note 
that we do not present the minimum values of the cross sections if it is found to be less than 
10 “'^ fb. Also, we only present single numbers for those cases where we have very few allowed 
points in the parameter space and no presentable variation in rates is observed for them. It is 
well known that the signals in the 77 , rr channels are comparatively clean due to the efficient 
tagging of photon and improved detection of r leptons either in jets or leptonic final states. 
For higher values of mjjj, say for 110 GeV, the dominant cross section is found to be again for 
the bb final state values ranging over from a few fb to a few thousands fb, whereas in the case 
of the TT final state, the maximum value lies within a few 100 fb. Interestingly, note that the 
rate of the Hi production decaying in the Hi —?• gg channel is reasonably large for all masses, 
unlike the situation in the SM as well as in the MSSM, where it is very much suppressed. 
Hence this channel can have the potential candidate to distinguish the NMSSM from other 
SUSY models, although elimination of QGD background is a non-trivial task. Remarkably, 
we observe that at higher values of mui for which mAi < mHi/2, Higgs to Higgs decays, such 
as Hi —)• AiAi open up with a maximum rate found to be 0(1000 fb) which is not shown 
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Table 2: Production rates(cross section x branching ratio)( in fb) of Hi at 13 TeV LHC 
energy for case(A). 


here. In this case, of course, the possibility of observing the signal of Hi depends on Ai decay 
pattern and the measurement can provide the information about triple Higgs bosons coupling. 
It is also observed that the rate in the WW channel at the higher masses range(~ 120 GeV) 
is non negligible, and can vary between the range of 0(1 - 100 fb). In addition, possibility of 
the Hi — 7 - XiXi (invisible decay) exists for higher mHi ( ^ 100 GeV) with a very tiny rate in 
a certain region of parameter space. 

Similarly in Table we present the minimum and the maximum rates of Ai in any given 
decay channel through both the production modes for a range of masses 22 - 249 GeV for 
the scenario represented by the case(A). The dominant rates of Ai also appear to be in the 
bb decay channel. On the other hand, rr, 77 channels predict a range of sizeable rates, in 
particular for mA^ ~ 100 GeV or less. For higher masses of Ai, the rates for rr, 66 channels 
are reasonable. Decays of Ai in various SUSY particles, in particular, for higher Ai masses 
also appear with a few fb cross section which are not presented here. In Table as before, 
we present the rates of the Ai production for a mass range from 17 — 249 GeV in multiple 
channels for case(B) where Hi is the SM-like Higgs. Note that as before the experimentally 
clean 77 and rr channels provide a sizeable rates for lower values of mAi, whereas for higher 
values of mAj, depending on the regions of the parameter space, the rates for rr, 66 channels 
appear to be feasible for exploring Ai. The cc, gg channels are not very promising due to tiny 
rates and presence of large QGD background, except for the lower masses where the rates 
are sizeable. As observed before for the Hi case, the bh channel has the most dominant rates 
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Table 3: Production rates(cross section x branching ratio) (fb) of Ai at 13 TeV LHC energy 
for case(A). 

for all masses of Ai. The production rate of heavier non-SM-like Higgs boson H 2 in various 
final states are presented in Table for the case(B). For the lower mass range, just above 
the SM-like Higgs boson mass, the dominant rate is due to the bb mode as expected followed 
by the same in rr channel which are not shown. For higher mass range around ~ 200 GeV, 
the rates for the channels such as, WW, AiAi and HiHi are quite large. It is also observed 
the rate in the invisible decay mode {H 2 —)• XiXi) is not negligible in a certain region of 
parameter space. Furthermore, for very high values (mj/j >200 GeV), Higgs to Higgs decay 
modes appear to be the dominant ones and provide a moderate rates ranging between a few 
fb to a few tens of a fb. 
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Rates 

17 

30 

59 

ruA^iGeV) 

100 

138 

181 

220 

249 

Rbbigg) 

36253 

10164 

.01 - 686 

6.3E-4 - 99 

4.5E-4 - 29 

7 

0.55 

4.8 

Rbiibb) 

202 

1058 

0.004 - 263 

0.001 - 123 

0.001 - 33 

28 

1.27 

1.5 

Rcc{gg) 

39 

6.4 

0.25 

0.45 

0.07 

0.022 

1.2E-3 

0.18 

RcAbb) 

0.21 

0.66 

0.1 

0.03 

0.008 

0.002 

0.0001 

0.0013 

Rrrigg) 

2421 

833 

.001 - 62 

10 

3.12 

0.81 

0.07 

0.58 

Rrribb) 

13.4 

87 

4E-4 - 24 

lE-4 - 13.6 

lE-4 - 3.7 

3.2 

.15 

.184 

Rn-Agg) 

8.8 

3.0 

0.22 

0.035 

.011 

2.8E-3 

2.5E-4 

2.1E-3 

R^p.{bb) 

0.05 

0.31 

0.1 

0.05 

0.013 

0.011 

5E-4 

6E-4 

-^77 { 99 ) 

0.3 

0.2 

4E-3 - 0.11 

0.042 

0.058 

0.017 

8.4E-3 

0.035 

Rj~^(bb) 

2E-3 

0.01 

2E-3 - 0.058 

lE-4 - 0.06 

0.04 

6E-3 

3E-4 

2.1E-3 

Rggigg) 

305 

46.3 

2.5E-3 - 1.7 

0.95 

0.22 

0.1 

6E-3 

0.97 

Rgg{bb) 

1.7 

4.8 

0.66 

0.15 

0.025 

0.014 

7E-4 

7E-3 


Table 4: Production rates(cross section x branching ratio) (fb) of Ai in various decay channels 
at 13 TeV LHC energy for case (B). 


Rates 



mH^iGeV) 



137 

200 

300 

400 

Rbbigg) 

32 

1.1 - 7.68 

2.5E-4 - 0.18 

1.6E-4 - 0.05 

Rb-bibb) 

78 

7.4E-3 - 1.5 

9.1E-4 - 0.066 

1.5E-3 

R^oxoigg) 

- 

30 - 58.32 

5.0 

7E-4 

R^o^o (bb) 

- 

0.19 - 11.51 

3.5E-4 - 1.8 

lE-4 

Rwwigg) 

0.1 

19 - 149 

3.1 

lE-3 - 13 

Rwwibh) 

0.24 

0.98 - 3.7 

5.6E-4 - 0.28 

2.2E-4 - 0.26 

Ra^Ai {99^ 

- 

602 

37 

0.35 - 24 

RAiAiibb) 

- 

4 

3 

4.5E-2 - 0.95 

RhiHi igg) 

0 

0 

2.8 

3.6E-4 - 7.5 

RhiHi (bb) 

- 

- 

0.4 

0.15 


Table 5: Production ratesfcross section x branching ratio)(fb) of H 2 at 13 TeV for case (B) 





5 Higgs signal 


The discovery potential of the non-SM-like Higgs boson depends on the underlying production 
mechanism along with its subsequent decay spectrum and the resulting rate. Already many 
studies exist in the literature presenting the detection prospect of finding the non SM Higgs 
bosons in the NMSSM exploiting, in particular popular bb, rr modes producing Higgs bosons 
directly or indirectly via Higgs to Higgs decays (for a review and details see Refs. [14,85,102 
and references therein). However, here we try to emphasize the detection prospect of the 
non-SM-like Higgs signal via enhanced rate in the 77 final state and as well as the other final 
states such as gg, cc. 

As we have noticed in the previous sections, the HR in 77 mode of Ai and Hi can be very 
large, in particular for Ai, it is closer to 90-100% in a certain region of parameter space in 
the scenario case(B), see FigjO [18,23,24, 103 


and similarly the BR of SM-like Higgs boson 


of mass 125 GeV in the di-photon is also enhanced in comparison to SM Higgs boson [104 
The BR of Hi —)■ 77 is also not negligible and is observed to be around 2-3% for mni masses 
~60 GeV or above as seen in the Figj^and the Table 2 [105 . For example, for ~ 68 GeV, 
(see Table [^, at 13 TeV LHC energy, one can have maximum about few thousands events in 
77 channel for the integrated luminosity 100/6“^ For the same scenario, in case of Ai, for 60 
GeV mass, one can expect about fewer event ~10-15 (see Table [^ for the same integrated 
luminosity. Here, the large BR of Ai into two photons channel does not yield a large rate 
because of suppressed production cross section. As is well known, looking for a signal in the 
di-photon final state is very promising as photon is a very clean object to tag experimentally. 
It is remarkable that neither in the SM nor in the MSSM, the decay rates of the Higgs 
bosons in the di-photon channel is comparable to what is in the NMSSM. For instance, in 
Refs. [^ 103j it is shown that for the Hi, the relative signal rates i.e production cross sections 
times BR(Hi —)• 77 ) are almost larger by almost a factor of 5-7 what is predicted by the SM 
for mill ~ lOOGeV and similar pattern is also observed even for miij > lOOGeV case 


106 


In 

this current study we focus on the di-photon channel for a very low mass region(<125 GeV) 
of non SM-like Higgs bosons which are accessible at the current LHC experiments. After 
performing a very systematic scanning of the parameter space we present the minimum and 
maximum rates of di-photon production for a given mass of non-SM-like Higgs. In addition 
we also discussed the variation of the BR of the di-photon channel for this mass range along 
with other decay channels. It is worth to have this information in order to estimate the 
sensitivity of non SM-like Higgs production for a given integrated luminosity. Evidently, this 
diphoton channel appears to be a smoking gun signal to look for light neutral non SM Higgs 
bosons. Although some simple parton-level analyses were already carried out in the past to 
explore the detection possibility in this di-photon channel 106 , |107 , more detail analysis are 
required in the present context. The estimation of the corresponding SM background in this 
di-photon channel is also a challenging task, in particular, for the low mass Higgs bosons. We 
postpone this analysis to a future work [108 . 

Furthermore, for Higgs decays to two gluons. Hi —gg with a branching ratio of about 
50% or more, a reasonable number of events in this channel are expected for the 100/6“^ 
integrated luminosity (see Fig. [^ and Table [^ , in particular for the mass range ~60-100 GeV, 
this channel yields a huge number of events ( 0 ( 1000 )) even for an 1 integrated luminosity as 
low as 10/6“^. Clearly, this gg final state provides an alternative option to detect the NMSSM 
light neutral Higgs bosons at the LHC by probing the invariant mass of a pair of jets, which 
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is not a easy task due to the presence of enormous QCD background. Perhaps, in the Hi+1 
jet final state, where an extra jet comes from radiation followed by Hi —)• gg decay, one can 
tag two jets to reconstruct the Higgs mass using the jet substructure method which is a very 
efficient technique 109 to probe boosted objects. Of course, the feasibility of this channel 


can be understood only after carrying out a detailed simulation including QCD background. 
Possibly, this gg channel provides a direct opportunity for a direct measurement of effective 
Higgs coupling to gluons directly [110 . In addition, the decays of the Higgs bosons to SUSY 
particles, such as, </>—?• if kienamtically allowed and have a reasonable HR, can also 

lead to interesting signals, particularly for leptonic decays of sparticles. Essentially, sparticle 
decay channels open up for values of niHi ~200 GeV. These decay channel may also 

contribute to the invisible width of the Higgs boson. Finally, Higgs to Higgs decay channels 
(Eq 3.26) give rise to many more detectable final states depending on the decays of Hi and 
Ai. A detail investigation of these channels and measurements of triple Higgs couplings may 
reveal the model structure of the NMSSM. In Table for the sake of illustration we present 
for few cases some of the relevant parameters where the rates are reasonably large for a given 
Higgs boson mass and decay channel. 



PI 

P2 

P3 

P4 

Case 

A 

A 

B 

B 

Mhi (GeV) 

13 (i=l) 

67 (i=l) 

200 (i=2) 

400 (i=2) 

Rxx(gg){fb) 

144927(66), 13838(rT) 

93 ( 77 ), 2881((75), 3737(cc) 

602(AiAi), U9{WW) 

24(AiAi), 13.03(IPTT) 

A, K 

0.363, 0,148 

0.325, 0.167 

0.544, 0.265 

0.64, 0.317 

tan P jj, 

9.27 ,165.45 

8.36 ,229.02 

2.45 ,192.7 

2.33, 401.30 


1440.7, -273.95 

1658.38,-450.35 

508.06, 3.69 

530.27,-35.35 

Mq3, Mus 

2828.25, 358.58 

1847.6, 2950.2 

2593, 1686.8 

1138.32, 1427.90 

At 

-2574.43 

-2779.78 

-3000.0 

-3016.5 

A/i, M 2 

89.18, 335.99 

75.9, 414.18 

92.89, 816.69 

80.6, 114.3 


P5 

P6 

P7 

P8 

Case 

A 

B 

A 

B 

Mai (GeV) 

22 

60 

100 

249 

Rxx{gg){fb) 

72336(66), 5939(rT) 

0 . 112 ( 77 ) 

226(66), 23(rr) 

0 . 035 ( 77 ) 

A, K 

0.368, 0.103 

0.328, 0.434 

0.677, 0.182 

0.499, 0.476 

tan P fi 

8.56 ,154.74 

5.61 ,128.7 

7.66 ,193.0 

1.55, 1859.0 

Ax , A^ 

1432.96,-0.486 

553.3, -1.53 

1464.15,-81.22 

872.53,-13.41 

Mq3, Mu3 

2237.9, 829.6 

2744.15, 1441.90 

1962.45, 2666.53 

1289.96, 1165.72 

At, Ab 

-2516.9 

-3914.36 

1850.3 

-1617 

Ml, M 2 

396.07, 319.15 

98.19,503.6 

410.8,160.43 

56.2, 574.26 


Table 6: Parameter sets(Pl - P8) for a given Higgs boson mass and decay channel for both 
the cases A and B corresponding to the rates Rxxij^ fb). These points P1-P8 are selected 
from Tabe 2-5. 
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6 Summary 


The recently discovered Higgs particle can be accommodated in the framework of the NMSSM 
easily without much adjusting the input parameters unlike in the case of the MSSM. Scanning 
the parameter space taking into account all constraints, we found the scenario where one of 
the Higgs bosons is SM-like and, some of the singlet like neutral Higgs bosons can be rather 
very light with suppressed couplings to fermions and the gauge bosons. In this work we 
focus our attention only to light neutral two CP even Higgs bosons and one CP odd Higgs 
boson out of which one the CP even Higgs bosons is SM-like. We discuss the masses and 
the couplings of the non-SM-like light neutral Higgs bosons and their BRs in various decay 
channels which predict the Higgs phenomenology at colliders. We present the BRs of the non- 
SM-like Higgs bosons in all accessible decay channels in reference to the allowed parameter 
space. We present the ranges of over all rates of production for a given non SM Higgs mass 
and for a given decay channel in the allowed region of parameter space, which demonstrate 
the potential final states that can be probed in the search for non SM signal at the LHC. We 
found some of the characteristic decay channels which are less dominant in the the SM and in 
the MSSM could play important roles in such searches. For example, in certain region of the 
parameter space, a large BR of Higgs decay to two photons is found to be quite interesting 
and can provide not only a characteristic robust signal of the NMSSM, but can also help to 
discriminate the NMSSM from the MSSM. Although this observation is also made in previous 
analyses, as discussed in sec. 5, but here we present the rates in this di-photon channel more 
systematic way for various lower masses of Higgs bosons corresponding to the center of mass 
energy at 13 TeV. These information may be very useful in searching for the non SM-like 
Higgs bosons at the LHC. Higgs decay to two gluons also turn out to be an interesting mode 
which one can exploit to study the the NMSSM Higgs sector, though care must be taken to 
deal with an enormous QCD background. Moreover, further studies of Higgs to Higgs decays 
via triple Higgs bosons couplings and probing the associated signals in various final states 
could uncover the underlying dynamics of the model. 
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